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Abstract. We demonstrate the role of the mass asymmetry on the balance energy (Et a i) by 
studying asymmetric reactions throughout the periodic table and over entire colliding geometry. 
Our results, which are almost independent of the system size and as well as of the colliding 
geometries indicate a sizeable effect of the asymmetry of the reaction on the balance energy. 



1. Introduction 

Reaction dynamics has always captured a central place in nuclear physics mainly due to the 
wider domain of physics it caters. Right from the low energy (where fusion-fission takes place) 
to the ultra-relativistic energies (where sub-nucleonic degrees of freedom are dominant) a large 
number of new phenomena has been observed/predicted. The reaction dynamics in heavy-ion 
collisions at intermediate energies has been used extensively during the last three decades to 
understand the nature of nuclear matter at extreme conditions of temperature and density and 
the nature of nuclear equation of state (EoS) [H [2], 13 HI E], 13 [7] • The phenomena that are mainly 
observed in this energy range are mainly multifragmentation, collective transverse flow, particle 
production etc. 

In the search of nuclear equation of state as well as of nuclear interactions and forces, collective 
flow has been found to be of immense importance [H [2j El HI El El [7J . Among collective flow, 
transverse in-plane flow enjoys special status. During the last two decades, much emphasis 
has been put on the study of collective flow [H El El El El El Ej- Lots of experiments have 
been performed and number of theoretical attempts have also been employed to explain and 
understand these observations. As reported by [HI E] for the first time and later on by many 
others, collective flow is negative at low incident energies whereas it is positive at a reasonable 
higher incident energies. At a particular incident energy, however, a transition occurs. This 
transition energy is also known as energy of vanishing flow or balance energy (-E& a j). This 
balance energy (Ef, a i) has been subjected to intensive theoretical calculations using variety of 
equations of state as well as nucleon-nucleon cross-sections [U [21 SI O El El El E] - This 
also includes the mass-dependence of E^i which have been reproduced successfully by various 
theoretical models [U El II El El E] • 

Interestingly, most of these studies take symmetric or nearly symmetric reactions into account. 
Recently, FOPI group studied the flow for the asymmetric reaction of 40 Ca+ 197 Au [lOUllj . They 
noted that the flow in the asymmetric collisions is a key observable for investigating the reaction 
dynamics. In contrast to symmetric collisions, where center of mass is one of the nucleus, this 



quantity is not known a priori in asymmetric nuclei experimentally. Later on, FOPI conducted 
experiment on 5S Ni and 208 Pb jlQ|, lllj . In other class of studies, the total mass of the system was 
kept fixed as 96 units whereas charge was varied [12]. Theoretically, recently Kaur and Kumar 
|13j conducted a complete study of the multifragmentation by varying the asymmetry of the 
colliding nuclei. Asymmetry parameter r\ is defined as (At — A p / At + Ap); where At and Ap 
are the masses of the target and projectile, respectively. All these attempts point towards a need 
for the systematic study of the disappearance of flow for asymmetric colliding nuclei. Further 
as noted, asymmetry of a reaction plays dramatic role in heavy-ion collisions [144 [To! [TB| I17| . 
This happens because excitation energy in symmetric colliding nuclei leads to larger compression 
while asymmetric reactions lack the compression since large part of the excitation energy is in 
the form of thermal energy. 

Note that some isolated studies with asymmetric nuclei are already done in the literature 
where the reactions of 20 Ne+ 12 C, 20 Ne+ 27 Al, 20 Ne+ 63 Cu, 58 Ni+ 12 C, 64 Zn+ 27 Al, 1 H+ 197 Au, 

12 C + 197 Au> 197 Au + 12 C; 40 Ar + 207 p h etc are taken into account [1^ [Til [JHl [IS IZOl [211 E21 [23l 

[24l 1251 [26] . Interestingly, none of the studies focus on the i?& a ; for asymmetric colliding nuclei. 
To address this, we here present a systematic study of the E^ai as a function of asymmetry of 
the colliding nuclei. While total mass of the reactions remain fixed, the asymmetry is varied 
by transferring the neutrons/protons from one nucleus to other. The QMD model used for the 
present analysis is explained briefly in the section 2. Results and discussion are presented in 
section 3 followed by summary in section 4. 



2. Description of the model 

The quantum molecular dynamics model (QMD) simulates the reaction on an event by event 
basis [2]. This is based on a molecular dynamics picture where nucleons interact via two and 
three-body interactions. The explicit two and three-body interactions preserves the fluctuations 
and correlations which are important for iV-body phenomenon such as multifragmentation [2]. 

In the QMD model, the (successfully) initialized nuclei are boosted towards each other with 
proper center-of mass velocity using relativistic kinematics. Here each nucleon a is represented 
by a Gaussian wave packet with a width of \f~L centered around the mean position r a (i) and 
mean momentum ^(t) [2]: 

Mr,P,t) = -^_ e [-{--W} 2 /4L] e [^ W -^]. (1) 
(2vrL) d/4 

The Wigner distribution of a system with At + Ap nucleons is given by 

f(r,p,t) = ^et-l-^Wi'/^eh^W} 22 ^ 2 ]^ (2) 

with L = 1.08 fm 2 . 

The center of each Gaussian (in the coordinate and momentum space) is chosen by the Monte 
Carlo procedure. The momentum of nucleons (in each nucleus) is chosen between zero and local 
Fermi momentum [= y/2m a V a (r); V a (r) is the potential energy of nucleon a ]. Naturally, one 
has to take care that the nuclei, thus generated, have right binding energy and proper root mean 
square radii. 

The centroid of each wave packet is propagated using the classical equation of motion [2] : 



where the Hamiltonian is given by 



a " 

Our total interaction potential V tot reads as [2J 

-ytot yLoc _|_ -yYuk _|_ ^/Co«Z _|_ -yMDI ^g^ 

with 

= t 1( 5(f Q - f^) + i 2 <5(r Q - r /3 )5(r a - f 7 ), (7) 

F y„fc = ^e-l^-r.l/m/ _ ?fj y m) > (g) 

with m = 1.5 fm and £3 = -6.66 MeV. 

The static (local) Skyrme interaction can further be parametrized as: 

tLoc _ n (P \ , a(P N " 



Here a, /3 and 7 are the parameters that define equation of state. The momentum dependent 
interaction is obtained by parameterizing the momentum dependence of the real part of the 
optical potential. The final form of the potential reads as 

U MDI « t A ln 2 [h{p a - ppf + l]5(r a - r p ). (10) 

Here = 1.57 MeV and % = 5 xlO - MeV~ 2 . A parameterized form of the local plus 
momentum dependent interaction (MDI) potential (at zero temperature) is given by 

U = a f £\ +p (P-Y + Sln 2 [e(p/ Po ) 2 / 3 + l]p/p . (11) 

The parameters a, /3, and 7 in above Eq. (11) must be readjusted in the presence of momentum 
dependent interactions so as to reproduce the ground state properties of the nuclear matter. 
The set of parameters corresponding to different equations of state can be found in Ref. |2j. 

3. Results and Discussion 

For the present study, we simulated various reactions for 1000-5000 events in the incident energy 
range between 90 and 500 MeV/nucleon. In particular, we simulated the reactions of § 7 0+ii Na 
(77 = 0.1), fN +f| Mg (77 = 0.3), fB +f A Si (7/ = 0.5), and \IA +f| S (77 = 0.7) for total mass 
(A TOT ) = 40, UAr+^Ca (77 = 0.1), %Si+%Cr (r, = 0.3), l° Ne+f s Ni (77 = 0.5), and fB+^Ge 
(77 = 0.7) for total mass {A TOT ) = 80, |§Ge+gjZr (7/ = 0.1), ^Fe+lfCd (7/ = 0.3), f Q Ca+\fTe 
(77 = 0.5), and HMg+\fCe (77 = 0.7) for total mass (A TOT ) = 160, and {fCd+lfBa (77 = 0.1), 
USr+lfDy (7/ = 0.3), f%Ni+]fW (77 = 0.5), and HAr+lfPb (77 = 0.7) for total mass (A TOT ) 
= 240. The present study is for peripheral collisions (i.e. b/b max = 0.5). Note that in some 
cases, slight variation can be seen for charges. The charges are chosen in a way so that colliding 
nuclei are stable nuclides. A momentum dependent soft equation of state with standard energy 
dependent cugnon cross-section (labeled as SMD) is used for the present calculations. 

The balance energy (Et, a i) is calculated using the directed transverse momentum < P!~ ir >, 
which is defined as: 

(P* ir ) = ^£sign{y(7)} p x (i), (12) 
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Figure 1. The time evolution of the directed transverse flow < P£ ir > as a function of time 
for different mass asymmetries. Columns 1 and 2 from left are for system mass {Atot) = 40, 
whereas columns 3 and 4 are for system mass (Atot) = 240. Here we display the results for 
peripheral collisions at different incident energies using SMD equation of state. 



where Y(i) and p x (i) are the rapidity distribution and transverse momentum of i particle, 
respectively. 

In Fig. 1, we display the time evolution of the directed transverse flow < P^ ir >, for various 
asymmetric reactions at different incident energies with system mass (Atot) = 40 and 240 
units. From the figure, it is clear that transverse flow is always negative during the initial phase 
of the reaction irrespective of incident energy, system mass, and asymmetry of the reaction. 
This shows that the interactions among nuclei are attractive during initial phase of the reaction. 
These interactions turn repulsive depending on the incident energy. The transverse in-plane flow 
in lighter colliding nuclei for all values of r\ saturates earlier compared to heavy colliding nuclei. 
A sharp transition from negative to positive flow is seen for lighter nuclei at lower asymmetry. 
For the heavier nuclei at all asymmetries; the transition is gradual. From the figure, one can see 
that nearly symmetric reactions (r/ = 0.1 and 0.3) respond strongly to the change with incident 
energy compared to highly asymmetric reactions. The cause behind is that with the increase 
in the asymmetry of a reaction lesser binary collisions take place resulting in lesser density and 
therefore, less response occurs for variation in the collective flow. On the other hand, for nearly 
symmetric reactions, the binary collisions increases linearly with incident energy, therefore, huge 
difference can be seen. 

It would be interesting to see how mass dependence character of the E^i behaves at a 
fixed asymmetry. In Fig. 2, we display the E^i verses combined system mass by keeping the 
asymmetry fixed. We notice that in all the cases, a perfect power law dependence (with power 
factor close to 1/3) can be seen for all asymmetries right from 0.1 to 0.7. The values of power 
factor are -0.34, -0.36, -0.35, and -0.34, respectively, for r\ = 0.1, 0.3, 0.5, and 0.7. Strikingly, the 
results are similar to the one as shown recently by Chugh and Puri |27| . There it was predicted 
that the role of MDI for symmetric systems is more dominant at peripheral geometries and the 
value of power factor is nearly equal to -0.35. From the figure, it is clear that all points are 




Figure 2. (Color Online) The Ej, a i as a function of system mass (Atot) of reacting partners for 
SMD EoS and for peripheral collisions. The results for different asymmetries r] = 0.1, 0.3, 0.5, 
and 0.7 are represented, respectively, by the half filled squares, circles, triangles and inverted 
triangles. Lines are power law fit oc A T . 
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Figure 3. The E^ a i as a function of system mass. Solid squares(circles) are for Soft lso (SMD) 
EoS. The top(bottom) left and right panels are for asymmetry parameter of 0.1(0.5) and 0.3(0.7) 
respectively. Lines are power law fit oc ATp OT . 



lying on the line, indicating that asymmetry plays a major role in the mass dependence. The 
variation in E^i for r/ varying between 0.1 and 0.7 is ~ 32, 35, 19, and 19 for Atot = 40, 80, 
160, and 240, respectively. This explains why in the earlier mass dependence studies [1], though 
average behavior was a power law, individual E^ai were quite far from the average law [HGiJIl]. 
There this happened because no control was made for the asymmetry of a colliding pair. Had 
all reactions been analyzed on a fixed asymmetry, one could have obtain perfect match with 
average power law. 

The effect of MDI on the E}, a i is shown in Fig. 3, where we display the E\, a i as a function 
of combined mass of the system with Soft* 50 (Soft 4SO represents the soft equation of state with 
isotropic standard energy dependent cugnon cross-section) and SMD EoS for r] = 0.1, 0.3, 0.5, 
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Figure 4. (Color Online) E}, a i as a function of asymmetry parameter rj is displayed in the upper 
panel, whereas the percentage difference AEb a i(%) as a function of system mass of reacting 
partners is shown in the lower panel. SMD EoS and an reduced impact parameter of 0.5 is used 
in both cases. The results for the systems having total mass Atot of 40, 80, 160, and 240 are 
represented, respectively, by the solid squares, circles, triangles and inverted triangles. Lines 
are to guide the eye. The results of the percentage difference for different asymmetries r/ = 0.1, 
0.3, 0.5, and 0.7 are represented, respectively, by the half filled squares, circles, triangles and 
inverted triangles. 



and 0.7. Lines represent the power law fitting and different symbols are explained in the caption 
of the figure. It is clear from the figure that as r\ increases, the suppression of E^ai for the lighter 
systems with MDI also increases while the increase in E^ a i for heavier nuclei remains almost 
same. Since transverse flow decreases with increase in asymmetry due to the decrease in nn 
collisions and Coulomb repulsions, therefore, the effect of the repulsive nature of MDI increases 
with increase in asymmetry mainly for the lighter masses. With the increase of asymmetry 
of the reaction, the role of MDI on the mass dependence of E^ a i is similar to its role for the 
symmetric reactions at semi-peripheral and peripheral geometries. 

In Fig. 4, we display E^ai as a function of 77 for a fixed mass equal to 40, 80, 160, and 240 
units (top panel). In agrement with all previous studies, E^ai decreases with increase in the mass 
of the system. This decrease has been attributed to the increasing role of Coulomb forces in 
heavier colliding nuclei. As discussed in earlier figures, a sizeable influence can be seen towards 
Ebai with variation of the asymmetry of a reaction. For lighter masses, the effect of the variation 
of 77 can result about 35 MeV change in the E\, a i- As noted in absolute terms, lighter nuclei are 
more affected compared to heavier ones. Overall, one sees that the effect of the asymmetry of a 
reaction is not at all negligible. It can have sizeable effect which goes as power law with power 
factor close to 1/3. We also display the percentage difference AEb a i(%) defined as AEbaiCVo) 
= ((£m°-£Cz°)/£Cz°)x100 as function of system mass (bottom panel). Very interestingly, we 
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Figure 5. (Color Online) The £y as a function of reduced impact parameter. The top(bottom) 
left and right panels are for Atot = 40(80) and 160(240) respectively. The results for different 
asymmetries r\ = 0.1, 0.3, 0.5, and 0.7 are represented, respectively, by the half filled squares, 
circles, triangles and inverted triangles. Lines are to guide the eye. 

see that the effect of the asymmetry variation is almost uniform throughout the periodic table. 
In other words, asymmetry of a reaction can play significant role in E^i and the deviation from 
the mean line can be eliminated if proper care is taken for the asymmetry of colliding nuclei. 

In Fig. 5, we display the Et> a i as a function of the impact parameter for different asymmetries. 
A well known trend i.e. increase in the E\, a i with impact parameter can be seen for all ranges 
of system mass. Further, as demonstrated by many authors, the impact parameter variation is 
less effected in the presence of MDI. The striking result is that the effect of mass asymmetry 
variation is almost independent of the impact parameter. 

4. Summary 

Using the quantum molecular dynamics model, we presented a detailed study of the balance 
energy with reference to mass asymmetry. Almost independent of the system mass as well 
as impact parameter, an uniform effect of the mass asymmetry can be seen at the energy of 
vanishing flow. We find that for large asymmetries, (77 = 0.7), the effect of asymmetry can 
be 15% with MDI and in the absence it can be 40%. This also explain the deviation in the 
individual E\, a i from the mean values as reported earlier. 
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